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Abstract. The system of N classical particles on the line with the Weierstrass 
p function as potential is known to be completely integrable. Recently D'Hoker 
and Phong found a beautiful parameterization by the polynomial of degree N of 
the space of Riemann surfaces associated with this system. In the trigonometric 
limit of the elliptic potential these Riemann surfaces degenerate into rational 
curves. The D'Hoker-Phong polynomial in the limit describes the intersection 
points of the rational curves. We found an explicit determinant representation of 
the polynomial in the trigonometric case. We consider applications of this result 
to the theory of Toeplitz determinants and to geometry of the spectral curves. 
We also prove our earlier conjecture on the asymptotic behavior of the ratio of 
two symplectic volumes when the number of particles tends to infinity. 



1. Introduction. Since its discovery in the middle of seventies by Moser and 
Calogero, the system of classical particles on the line with the Hamiltonian 



c3 ; N ? 2 2 N 



H = ^2 y — p ( qn ~ 9r )' where a 2 = ±1, m > 0, 



n,r=l 



has attracted a lot of attention from mathematicians and physicists. This sys- 
tem is completely integrable in the sense of Liouwille, i.e., it has N commuting 
integrals of motion. Algebro-geometrical methods in the study of this system 
were introduced by Krichever, [Kl]. He showed that the Riemann surfaces - 
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spectral curves, associated with this system are iV-sheeted covers of the elliptic 
curve of the potential, the Weierstrass p function. The geometry of such covers 
is fairly complicated and was not understood well until now. 

In the papers of Seiberg-Witten and Donagi-Witten, [SW,DW], this system 
was considered in connection with Super-Symmetric Yang-Mills quantum field 
theory. Degeneration of the elliptic potential into a trigonometric one is inter- 
preted by physicists as a classical limit of the quantum system. In such limit 
genus- N spectral curves degenerate into genus-0 curves. Recently D'Hoker and 
Phong, [DP], found a beautiful parameterization by a polynomial Q(p) of the 
space of spectral curves for particles with elliptic potential. In the trigonometric 
limit the zeros k n of the polynomial Q(p) = m N Y\n=i(P ~ ^n) describe in- 
tersection points of the limiting rational curve. We have found a determinant 
representation of the polynomial Q(p) in this limiting trigonometric case. Our 
description becomes completely explicit when the mechanical system is at the 
ground state: all particles at rest, placed equidistantly on the circle. In this case 
Q(p) can be computed explicitly using the determinant representation. This 
provides us with complete information about the spectral curve. In all other 
cases, such as an elliptic or a system not at the ground state, the curve can be 
viewed deformation of this simple case. 

In this paper we present a simple proof of the D'Hoker-Phong theorem, our 
determinant representation of Q(p) and some of its' various applications. The 
paper is organized as follows. In section 2 we present some well-known facts 
about particles with elliptic potential. In section 3 we give a simple proof of 
the D'Hoker-Phong theorem and prove our main result: an explicit formula 
for the polynomial Q(p). In the three subsequent sections we consider various 
applications of our result: to the theory of Toeplitz determinants (section 4), to 
the geometry of spectral curves (section 5) and finally to the asymptotic behavior 
of the ratio of two symplectic volumes, when the number of particles tends to 
infinity (section 6). 

2. Commutator formalism. Riemann surface. The iV-particle Hamilton- 
ian 

N 2 2 2 N 

n=l n,r=l 

produces the equations of motion 
• OH 

(2.1) q n = - — =p n , n=l,...,N, 

OPn 



p n = -—— =m a y^p'{q n -qr), n = l,...,N. 
dq n ^ 

For attractive particles the parameter a = 1 and for repulsive particles a = i; 
the parameter m > plays the role of mass. 
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For m = 2 the iV-particle system can be embedded into elliptic solutions of the 
Kadomtzev-Petviashvilli (KP) equation. The key step is the following theorem 

Theorem 1. [Kl]. The equations 



N 



n=l 
n 

ad y - + p(x - q n (y)) 



</> = 0, 



n=l 



have solutions of the form 



N 



ip(x, y, k,z) = J2 a n(y, k, z)Q(x - q n , z)e kx+a ' k y , 

n=l 
N 

i(j + (x, y, k,z) = J2 a n(y, k, z)Q(-x + q n , z)e~ kx - a ~ lk2y , 



n=l 



where* <fr(q,z) = ^z)a(q) e< '^ 9 > ^ an< ^ on ^/ ^ ^ e f unc ^ ons Qn(y) satisfy the 
system of equations (2.1). 

The operator L 2 with the elliptic potential 



N 



L 2 = d 2 x - u(x, y) = dl - p(x - q n (y)) 

n=l 

enters into the commutator formula 

[ad y -L 2 ,d t - L 3 ] = 0, 
with L 3 = d x — ^ud x — w. These imply that u satisfies the KP equation 

— o~ Uyy {Uf + ~uu x ~~~^a xxx ) x . 

For arbitrary mass m > 0, the system (2.1) is equivalent to the matrix equa- 
tion aL = [L, M], with N x N matrices M and L having entries 

L nr = o-p n 5 nr + m$(q n - q r ,z)(l - S nr ), 

M nr = \ -p(z) +m^2p(q n -q s ) \ 5 nr + m&(q n - q r , z)(l - 8 nr ). 



*a{z) denotes the Weierstrass function. See section 7 for definitions. 
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The matrix L can be simplified using a gauge transformation 

L = GLG~\ 

where G nr = e^ z)qn 5 nr and L nr = ap nr 5 nr +m$o(q n -q m , z)(l-S nm ), $ (q, z) = 
a(z)a(q) • ^ ne s P ec trum of L is preserved and this determines the curve 

T = {(k,z): R(k, z) = det(L + mk) = det(L + mk) = 0}. 

We denote by lower and upper Greek letters the points (k, z) on the curve T, for 
example, 7,11, etc.. 

Remark. An exchange of positions and velocities of particles labeled, say, n 
and r does not affect the curve. It simply permutes the n-th and r-th columns 
and rows of the matrix L. 

The following facts about the curve can be found in [Kl]. 
(i.) The determinant R(k, z) can be written in the form 

N 

R(k,z) = J2r n (z)k n , 

where r n (z) are elliptic functions of z. This means that V is an iV-sheeted 
covering of the elliptic curve To- 
(ii.) In the vicinity of zero 



R(k, z) = m 



n ■ (^+M o, + ...)]n[*-(4 + *«°' + ...) 



The points II n , n = 1, . . . , N above z = are called infinities. The infinity 111 
corresponds to the "upper sheet", where k(z) = + 0(1). 
(Hi.) In the elliptic case, for the generic configuration of iV particles the genus 
of the curve V is N. In the trigonometric case the curve is rational. 

The following lemma describes symmetries of the matrix L. We assume that 
2uj is real and 2u' is pure imaginary. 

Lemma 2. (i). a = 1. Let T\ is defined as: T\(k,z) = (k,z). Then 



(L + mk) (717) = (L + mk) (7). 
(ii). a = i. Let r» is defined as: Ti(k,z) = (—k, —z). Then 



(L + mk) (ra) = —(L + mkf (7). 
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Proof, (i). The statment follows from the identities 

(L + mk) nn (T 1 'j) =p n + mfc(7) = p n + mk{^) = (L + mk) nn {^), 

(L + mfc) nr (ri7) = m$(<? n - qf r ,z) = m$(q n - q r ,z) = (L + mk) nr (7). 

(ii). The proof is similar to that in the attractive case. □ 

Lemma 2 implies that in the attractive/repulsive case there exists an antiholo- 
morphic involution ti/tj on the curve V. 

For some special configurations of particles there is also an additional sym- 
metry. 

Lemma 3. Let r_ is defined as: T-(k,z) = (—k,—z). If p n = for all n = 
1, . . . , iV; then 

(L + mk) (r_7) = — (L + mk) T (7). 

Proof. It is enough to note that 
(L + m/c) nn (r_7) = -mk^) = — (L + mk) nn {^), 

(L + m/c) nr (r_7) = m$(q n - q r , -z) = -m$(q r - q n , z) = — (L + mk) rn {^). 
□ 

Lemma 3 implies that if p = 0, then there is an involution r_ on the curve. 

3. D'Hoker-Phong parametrisation. Determinant Formula. Following 
D'Hoker-Phong, [DP], let us introduce 

The function h\(z) = ((z) — has periodicity properties 

h\{z + 2uj) = h\(z), 

h 1 (z + 2u') = hi_{z) . 

Let us also introduce a multivalued function ^(7) on the curve Y as 

pin) = Hi) + h i( z (l))- 

The function ^(7) has a simple pole at LTi only and does not change under the 
shift z — > z + 2uj. It is defined up to an integer multiple of ^ which corresponds 
to the increment of h\ under the vertical shift z — > z + 2a;'. It is called quasimo- 
mentum, because ip(x, y, 7) defined in Theorem 1 has the Bloch property in the 
x-variable: 

ip(x + 2a;, y, 7) = e p{l)2uj 'ip (x , y, 7). 

So far we represented Y as a curve in the direct product of the k and z variables. 
Now we want to write an equation for Y in the direct product of p and z. 
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Theorem 4. [DP]. There exists a polynomial Q(p) of degree N, 

N N 



Qip) = m N J2p N ~ n Qn = m N H(p- k n ), 



n=0 n=l 

such that 

' z-d\ °° 



n I z-d \ oo 7 / \ 

Mid n! 



Proof. * Performing the change of variables £(7) = p(7) — /ii(z(7)) we have 
R(p — hi(z),z) = m N rin=i (P ~ Pn( z ))- The roots p n (z) do not change under 
the shift z — > z + 2cu. They transform as 

p Sr (z + 2a;') = p r (» - — 

under the vertical shift z — > 2 + 2o/; (si, . . . , sjv) is some permutation of sheets. 
Similarly the operator 

Mjg) 

does not change under the horizontal shift and for the vertical shifts we have 



Because of these, the polynomial 

N (—) 

Q(p, z) = m N J2 P N ~ n Qn(z) = ^MvR(p ~ h x {z\ z) 

n=0 Ul \~) 

has elliptic coefficients Q n (z). The coefficients of the polynomial on the right 
hand side 



R(p — hi(z), z) = m 



N 



have singularities at z = not higher than a simple pole. The operator 

^ (l^f ) /^(^j) does no ^ decrease the order of a singularity. Therefore, the 



*This short proof was communicated to me by I. Krichever 
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elliptic coefficients Q n (z) have a singularity at z = not higher than a simple 
pole, and these singularities do not depend on z at all. □ 

When we want to emphasize the dependence of all functions on modular 
parameter r we write h(z) = h(z\r), etc. In the trigonometric limit r — > +zoo 
an elliptic curve To degenerates into a sphere with two punctures, parametized 
by w = e~^. North ((3) and south (y) punctures correspond to w = oo and 
w = 0. The function ^(7) can be used as a global parameter on the limiting 
rational curve. The equation of the curve 



0i 



d 



2U! 



r Q(p\r) = 0, 



using the definition of 9\ can be written as 



Y {-l) n h n ^w~ n Q (p+ — - — r) 

\ 2uj uj I 



T = 0. 



In the limit r = +zoo it becomes 



/ ill \ _i / Z7T \ 



For each point v n e T, n = 1, . . . , N above w = there exists a point /? Sri G T, 
above w = 00 such that 



(3.1) 



p(u n )-- = k n =p((3 Sn ) + -. 



It is easy to check that v n and j3 Sn correspond to a simple crossing. 
Let us introduce the N x N matrix C such that 



£ nr = ap n 6 nr + mR(q n -q r )(l-6 nr ), where R(q) 



n 1 
2^imlgr- 



Theorem 5. The polynomial Q(p\r) for t = +ioo is given by the formula 

Q(p\ + ioo) = det (C + mp) . 



Proof. It is enough to prove that 



(3.1) lim det [L(z | r) + mp) 



= det (C + mp) , 
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lim h\{z\r) 



Z. 

- — 2 



(3.2) 

(3.3) lim 
After that, the D'Hoker-Phong formula taken at 2; 



Z7T 

2c? 



— (9 

= e 2 ^ . 



Q(p|r) = det ( £(z|r) + m(p — /ii (z) 



implies the statement, when r — > +zoo. 
In order to prove (3.1), note that 



and 



$ {q,z\T)\ z =T = 



lim 



1 



sn(g|r) 



-v 1 



r^+ioo sn (q\r) 



R(q). 



To prove (3.2), note that 
Therefore, 



^1 



when t — > +zoo. Finally, 





n 












in 


(- 


■) 


= i 




\2u 




2 =§ 


~2u_ 



nr 
2u 



(2n - 1) 



e ^^(2n-l)_ 



e 4 +0 



( e 4 )' 





Z7T 


*=5 


~2u;_ 



lim /ifc(z|r) 

r— >+ioo 

This completes the proof of (3.2). 

In order to prove the last identity (3.3) we use the previous result 



-a 



2lo 



T 



lim 



lim y 



h k {z\r) 
k\ 



(-d)> 



fc=0 v 



ITT Q 

e 2 ^ . 



2 fe=o 

In fact, only the first N + 1 terms in all sums matter, since the degree of each 
polynomial in p does not exceed N. This justifies interchange of the limit and 
sum. The last identity (3.3) is proved. □ 

There are two essentially distinct cases. For a = i the matrix C is skew-adjoint 
and all k n , roots of the polynomial Q(p), are pure imaginary. For a = 1 the skew- 
adjointness is lost, but similar to Lemma 3, if p = the skew-adjointness of £ is 
restored and all k n are pure imaginary. 
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4. Spectrum of the Toeplitz operator. For a finite number of repulsive 
particles (r < +zoo) there exists a point or set of points in the phase space 
where the Hamiltonian H(q,p), bounded from below, achieves its minimum. 
Such a point p n ,q n , n = 1, . . . , N is called a "ground state". In our case the 
Hamiltonian is rotationaly invariant and ground state is a one-parameter family: 
Pn = 0» q n - q n -i = if' for all n = 1, . . . , iV. 

Obviously, the matrix C at the ground state is Toeplitz and skew-symmetric. 
Let us intoduce e T = (e 1 , . . . , e N ) T , an eigenvector of C : (C + km) e = 
corresponding to the eigenvalue k and normalized such that e 1 = 1. 

Lemma 6. Let a = i. If the system is at the ground state, then the eigenvalues 
and eigenvectors of the matrix C are given by the formulas 

k n = JL(N + l-2n), n=l,...,N; 
loot 

and 

4 = e^0-i)(i-2n) ? j = l,...,N; n = l,...,N. 

To prove the Lemma we introduce some formulas from Fourier analysis on the 
additive group Gn = Z/NZ of the ring of residues modulo N. 

i. Characters \i G Gn are all defined as /j, n = e lf3n , (3 n = ^(n — 1), n = 1, . . . , N. 
The pairing is: for any s G Gn, H G Gn '■ (s, fx) = fi s . 

ii. For any complex function a = {a s , s G Gn} its Fourier transform is a(fi) = 

« s (s, A*), and the inverse transform is a s = Y <2(a0(~ s 7 A*)- 

seGjv n€G N 

Hi. For any complex functions a, b the Plancherel identity holds: Y a s~b s = 

sEG N 

7? a{y)b{n). 

Proof. To simplify calculations let m = ^ and also re-scale the spectral pa- 
rameter k — > 2^1^. General case easily follows from this. 

Step 1. Let e\ = (e 1 , . . . , e N ) T is an eigenvector corresponding to the eigen- 
value k. We will prove that e{\_ = (e N , . . . , e 1 ) 2 " is an eigenvector corresponding 
to the eigenvalue —k. Indeed, for any j 

^ C jr e r + + Cj :N -j+ie+~ j+1 + ke\ = 0. 

r^N-j + 1 



Then, 



^2 Cj,N-r+ie+ r+1 + C j}N -j+ie+ 3+1 + ke\ = 0. 



r^N-j + 1 
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At the ground state CN-j+i,N-r+i = —£j r - Introducing e™ = e^ _n+1 , we 
obtain 

^ CN-j+i, r e r _ - C N - j+1 je J _ - (-k)e^~ J+1 = 0. 

r^N-j + 1 

The statement is proved. 
Step 2. For j ^ r we have 

( A /9;\2, 1/2.1/2 

£ jr = mRty - q r ) = . 

Let D, C be the N x N matrices with entries 



D jr = V2i/j,y 2 8j r , 

k l 1 

Cj r = — — [Mj Sj r + _ (1 — fijr)- 
— ■ LL o LL f 

Then, £ + mk = DCD and the eigenvalue problem can be written as 

(C + mk)e = DCDe = Ce = 0, 

where e = De. For any j 

Introducing Uj = lij 1 *? and using jij/ ji r = ^j-r+i, we obtain 

+ V ^ = 0. 

This is a convolution equation which can be solved using Fourier analysis. In- 
troducing 

U{fJ>)= ^ u s (s,fj,), M(fj,)= (*>/*)» 

s€G n s6Gjv _{AT} ^ s + 1 

we arrive at the main equation 

U(fj)(2M(fj) - k) = 0. 
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Step 3. Note \x n = ^2 1 • For n > 2 we have 



MQMn) = J] 



s 



sea N \{N} ^ s+1 1 

^> _ 1 + 1 

Z / M s — 1 

s6G JV \{A r } P2 

n-2 

= m(pi)+ E E(^) fc 

seGjv\{iV} fc=o 

n-2 

= mm + E E + E 1S - 

fc=l seGiv\{Af} s6Gjv\{Af} 

If r N = 1 and r ^ 1, then r + r 2 + . . . r^ -1 = — 1. Therefore 

M(aO = M(a»i) + N + 1 - n. 

From this formula we see that the function M(p) takes different values at the 
different points of the character group Gn. Let k n = 2M( / u n ), n = 1, . . . , N 
in the main equation of Step 2. Then, k n = k± + 2(N + 1 — n), for n > 2. The 
eigenvalue k\ is the smallest and ki is the largest. Using Step 1: —k\ = k%, we 
finally obtain 

h = -N + l, 

k n = N + 3 - 2ra, w = 2, . . . , N. 

Step 4- Let, for any n 

Inverting the Fourier transform 



AteGjv 



Using the identities -u J = /u - 1 e J , e- 7 = ! e J ', we arrive at the formula 

e n = e^^ _1 ^ 3 ~ 2n ). Shifting n — >■ n + 1, completes the proof. □ 
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When the number of particles with m = 1 tends to infinity for the ground 
state with v = %^ = 1 the matrix C has the limit 




The corresponding Toeplitz operator has the symbol <p(t) = in — it, < t < 
2tv. To emphasize the dependence of k n on N we write k Uj N- For any continuous 
function / on iR 1 we have 

n=l 

Such distribution of eigenvalues is called canonical, see [W] . It is usually proved 
with the aid of Szego theorem. Our case is not covered by the conditions of the 
classical theorem, but the conclusion for this symbol is true also. 

5. Real ovals of anti-involution. Images of infinities. In the trigonometric 
case the curve Y is rational and p(Q), Q G Y can be used as a uniformization 
parameter. On the curve Y anti-involution ri/rj can be expressed in terms of 
the parameter p. 

Lemma 7. (i). a = 1. Then p{r\Q) =p(Q). 
(ii). a = i. Then p{tiQ) = —p(Q). 

Proof, (i). The statement follows from identity 



p(r 1 Q) = k(Q) + h 1 (z(Q))=p(Q). 

(ii). Similarly 

p(nQ) = -HQ) + hi(-4Qj) = -pW)- 

□ 

In the attractive corollary of the Lemma, we obtain, that = is a 

fixed oval of anti-involution t±, and in the repulsive case = is a fixed oval 
of anti-involution Tj. 
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Consider the repulsive case and let us look at the images of infinities IT. For 
n>2we have 

p(U n ) = U7) +M7)l 7 =n B = + ki n ] + 0(z) + -U=o = 4 0) - 

As is proved in Lemma 5 of [V2], all kn^ are pure imaginary and therefore all 
images of infinities LT2, . . . , II at lie on the fixed oval ?R,p = of anti-involution 7^. 
For n = 1 and 7 in the vicinity of LTi 

p( 7 ) = m + M7) = ^ + 0(1) + - = ^ + 0(1). 

z z z 

Therefore LTi corresponds infinity of the uniformization parameter p. 

For example, for iV = 2 at the ground state = and due to Lemma 6 we 
have the fig 1. 




not a ground state 

m 




FIG. 1 



FIG. 2 



When one pumps energy into the system the loops are perturbed, but kffl are 
still pure imaginary (fig 2). When the modular parameter r changes from +zoo 
(trigonometric case) to some finite imaginary number (elliptic case), the loops 
"move inside" and become fixed ovals around holes, as is shown in fig 3, 4. 
For N > 2 the picture is similar, but iV loops are present. 

The case of attractive particles is more complicated, because skew adjointness 
of the matrix C is lost. It is restored for p = 0. 
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ground state 



not a ground state 




T<1°° genus=l 



X< loo 



genus=2 



fig. 3 



fig. 4 



In the trigonometric case, for both repulsive and attractive particles the im- 
ages of infinities II2, . . . , IT at are moduli of the spectral curves. We will prove this 
statement in the elliptic case. First, we derive the trace formula relating total 
momentum P and kn\ n = 2, . . . , N. 

Computing ^ n res/c 2 cte = 0, we have J2n=i^ k n = 0- Using k n = 
— 1, n > 2 and k[~ r ^ = N — 1, we obtain Nk^ = J2n=i • On ^ ne other hand 
tr L = aP = —m J2n=i • By comparing the last two expression we arrive at 
the trace formula 



For another trace formula relating H and kn see [V2] . Now, we are ready to 
prove 

Lemma 8. In the elliptic case J n , n = 2,...,N, supplemented by Qn, the 
constant term of the polynomial Q(p), are moduli of the spectral curves. 

Proof. Using the representation of R(k,z) in the vicinity of zero and h±(z) = 
- + O(l) we have 



N 



mN 




where J„ = — k^ 



a(N-l)' 



71 = 2 




N 



z 
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FIG. 5 




Therefore 

N 

res R(p-h 1 (z),z) = -Nm N JJ (p - k^ ] ) . 

From the D'Hoker-Phong formula 

00 h ( \ N 

*n ( -* )fce(p) = - NmN ncp- 4 o) )- 

fc=0 ' n=2 

Since res/ifc(0) = for k even and res/ii(0) = 1 we can recover from J's the 
derivative of Q(p). The constant term Qn allows us to reconstract the whole 
polynomial Q(p)- □ 

6. Symplectic volumes. Let us intoduce in a general elliptic case the new 
Hamiltonian 

N a _2^2 N 



n=l n,r=l 



where p = p + — . In the trigonometric limit p(q) = R 2 (q)- 

Consider now the trigonometric case. Singularities of the curve V are simple 
crossings, as it is shown in fig 5. 

Form the identity (3.1) is follows that k{v n ) = k{(3 8n ) = k n . We will derive a set 
of new trace formulas. 

Lemma 9. In the trigonometric case the following formulas hold 



N 

P = I n , where I n = k n , 

1 ° 

n=l 

N 2 2 

H = t'm Where T 'n = -^~ k n- 

n=l 
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FIG. 6 



fixed ovals 




Proof. For Q(p) we have 



N 



Q(p)=m N l[(p-k n ) 

n=l 

= m N p N +p N - 1 [J2-^)+P N - 2 [ E *»*»') + 



K n=l 



, n,n' = l 



On the other hand we can write the polynomial Q(p) as a Fredholm expansion 
for the determinant given by Theorem 5 

Q(p) = det(£ + mp) = (mp) N + (mp)^" 1 ^ ap^j 

\ 

E a2 PnPn> ~ m 2 R(q n - q nl )R{q n > - q n ) + .... 



+ (mp) 



N-2 



. n,n' = l 



Comparing coefficients after simple algebra, we obtain the result. □ 

In the elliptic simple crossing becomes a handle, as shown in fig 6. 

In this case one can relate P, H to the periods of kdz, k 2 dz over a-circles of V. 

Lemma 10. In the elliptic case these formulas hold: 

N 



n=l 
N 



where I r , 



-7 / k[z) ^ 



2u' 



,2^2 



n=l 



where lL = ^lk 2 (z)^. 
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Proof. * Let c n be small contours around Il n , n = l,...,iV. Using Cauchy' 
theorem 

N „ AT 

^ / fc(z)fei(z)dz+5^ / A;(2)(/i 1 (2 + 2w / )-/ii(^))^ = 

n=l/ n=l/ 

Note that h\{z + 2a/) — h±(z) = — ^ and the second sum reduces to 

N . 

n=1 a 

To evaluate the first sum we use the asymptotics of k near infinities IT: 

J k(z)h 1 (z)dz = J + k ( V + 0(z) \ x Q + 0(zfj dz = 2mk$\ 



As a result we derive 

JV N „ , 

2a;' 



$>i 0) = £ 

n=l n=l„ 



On another hand we have 

TV 



trL = aP = -m^k^. 



n=l 



The first formula is proved. 

To prove the second trace formula one should integrate k 2 dz over the same 
contour. □ 



The meromorphic 2-forms 

JV N 

a 



tt = ^ dk (ln) A dz(^ n ) = ^2dI n A d(f) n 

n=l n=l 
2 2 N N 

^ = ^ ^ ^ = J2dI' n A dtp 



n=l n=l 



on the space of spectral curves can be written following [KP,K2] in a different 
way, namely 



*For another proof see [DP]. 
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Theorem 11. In the elliptic case the meromophic 2-forms can be written as 

N 



O = res n n < a*(5L — mk) A 5a > dz, 



n=l 
2 9 N 



Q' = m ° resu n k < a*(5L — mk) A 5a > dz. 



2 



The form O can be written as dq n A dp n , but for the higher form O' an 
explicit expression is not known. 

In the paper [VI] we considered a Jacobian between symplectic volumes 

N N 

An yj dl n 

Jac = — = || — . 

AO' n=l 01 * 

The quantity A^ -1 log Jac , when iV — > oo, has a meaning of entropy or, in prob- 
abilistic terms, rate function, see [L,Va]. 

Lemma 12. For repulsive particles in the trigonometric case at the ground state 
the following asymptotics hold * 

1 2 

— log Jac = log N + \ogv + log h o(l), when N — > oo; 

AT 7rm 

where v = 

Proof. Consider the case when A" is divisible by 4. At the ground state Lemma 
6 implies 

Therefore, 

N dl f 2 \ N N N 
JaC = naif = VW n AT + l_2n 

and 

1 2 2 

— log Jac = logw + log h log N — — log(A" — 1)!! 

N ' rmv N 

Using the formula 

— log(AT- 1)!! = log AT - l + o(l) 



*This formula was first conjectured in [VI]. 
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we have 

1 2 
— log Jac = logv + log h 1 + o(l). 

N run 

The remark of section 2 implies that in fact we consider action-angle variables on 
the fundamental domain q\ < q 2 < . . . < qN of the configuration space [0, 2u) N . 
Therefore the Jacobian has to be multiplied by N\. This and N^logNl = 
log N — 1 + o(l) imply the statement. 

The case of general N can be considered similarly. □ 

7. Appendix. The Weierstrass function a{z) has periods 2uj and 2o/ and is 
defined as 

„ w =zn{(i-5)«P(i+^)}, 

where u> = 2un + 2u'n' and Yl is taken with n, n' e Z 1 ; n 2 +n' 2 > 0; a(z) has 
zeros at the points of the lattice 2um + 2uj'n' . In the vicinity of zero it has the 
expansion 

a(z) =z + 0(z 5 ). 
Under the shift z — > z + 2u the function a(z) transforms as 

a(z + 2uj) = —cr(z) ex.p(2i](z + a;)). 

For the shift z — > 2: + 2a/ one should replace 2?y by 2?/ in the formula above. The 
functions ((z) and p(z) are defined as: 

= T~ lo § P(*0 = ~T"C(^)- 

They have respective expansions at zero of the form 

C(*) = - + 0(A p(z) = ^ + 0(^)5 
and periodicity properties 

p(^ + 2o;) = C(* + 2a;) = C(z) + 2 V . 

For the shift z — > ^ + 2a;' one should replace 2?y by 2?/ in the formula above. The 
Legendre's identity is r/o/ — r/'o; = 

The Jacobi functions 6> , #i are defined as 

n= — oo 
+ oo 
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where h = e inT . The relation between #1(^7) and a(z) is: 

a(z) = 777 ; 6>i — e 2 ». 
v ; ^(0) V2w/ 

The Jacobi sn (2) is defined as 

For more information see [HC]. 
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